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METHOD AND APPARATUS FOR IMAGING INTERNAL STRUCTURES OF 
TRANSPARENT AND TRANSLUCENT MATERIALS 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This application claims the benefit of the filing of U.S. Provisional Patent Application 
Serial No. 60/458.111. entitled "Method and Apparatus for Imaging Internal Structures of 
Transparent and Translucent Materials", filed on March 26. 2003. and the specification thereof 
is incorporated herein by reference. 

STATEMENT REGARDIN G FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT 
The U.S. Government has a paid-up license in this invention and the right in limited 
circumstances to require the patent owner to license others on reasonable terms as provided 
for by the terms of Contract No. DMI-0214911 awarded by the U.S. National Science 
Foundation (NSF). 

BACKGROUND OF THE INVENTION 
Field of the Invention (Technical Field^ : 

The present invention is of differential spectral interferometry (DSI). a novel method for 
biomedical and materials imaging, which combines the high dynamic range of optical 
coherence tomography (OCT) with Inherently parallel low-bandwidth image acquisition of 
spectral interferometry (SI). DSI efficiently removes the deleterious DC background inherent in 
SI measurements while maintaining the parallel nature of SI. The invention is demonstrated for 
both synthetic and biological samples. Because DSI is mechanically simpler than OCT. while 
preserving the low bandwidth, parallel nature of SI. it is competitive with OCT for biomedical 
applications in terms of image quality and acquisition rate. 



Background Art : 

Note that the following discussion refers to a number of publications by author{s) and year 
of publication, and that due to recent publication dates certain publications are not to be 
considered as prior art vis-a-vis the present invention. Discussion of such publications herein is 
given for more complete background and is not to be construed as an admission that such 
publications are prior art for patentability determination purposes. 

Optical coherence tomography (OCT) and spectral interferometry (SI) are low coherence 
optical imaging methods related to each other through the Fourier transfomi (FT). OCT. which 
takes measurements in the time domain, provides good dynamic range and image resolution. 
It has been demonstrated that OCT can be successfully used for sub-surface biomedical 
imaging at nearly video rate. A drawback of OCT is the requirement for a high-speed optical 
delay line in the reference arm. SI makes measurements in the frequency domain, taking 
advantage of inherently parallel spectral data acquisition and an optical arrangement without 
moving parts or high-speed modulation. Recently. Applicants reported on successful use of SI 
for 3D-imaging in Xenopus tadpoles. A.B. Vakhtin, et al., "Differential spectral interferometry: 
an imaging technique for biomedical applications". Opt. Lett. 28. 1 332 (August 2003); and A. 8. 
Vakhtin. et al. "Common path interferometer for frequency-domain optical coherence 
tomography," AppL Opt. 42. 6953 (December 2003). 



The main drawback of SI is the DC character of the measurements, which results in a 
relatively low imaging dynamic range. A strong background associated with the spectrum of 
the light source itself (DC term) and the interfering waves scattered from different surfaces 
within the sample (autocorrelation terms) are inherently present in SI images. It is the DC and 
autocorrelation terms that keep the dynamic range of SI far below the theoretical shot-noise 
limit, thus limiting the usefulness of SI for biomedical imaging. To eliminate these two artifacts 
of SI. M. Wojtkowski. et al.. "Complex spectral OCT in eye imaging". Opt. Lett. 27. 1415 
(August 2002). recently suggested an approach that involves complex five-frame phase and 
amplitude reconstruction. The method improves the signal-to-noise ratio and doubles the 



depth range by using both positive and negative optical path differences for imaging; however, 
it requires exceptional stability of the object — mO within the complex-method data 
measurement time scale, which Is 6-120 s. 

A major advantage of the high-speed differential detection of the invention is that it 
rejects 1/f noise. In the case of 1/f noise, the magnitude of low frequency components Is much 
larger than the magnitude of high frequency components. Thus, when the detection bandwidth 
allows low frequency components to pass, the noise levels are higher simply because the 
amplitude of the low frequency components is higher. 

Therefore, it Is imperative to reject low-frequency components, which Is not taught by 
either U.S. Patent No. 5.565.986. to Knuttel or U.S. Patent No. 6.377.349. to Fercher. While 
they do demonstrate that a dither In the reference arm can be added to enhance data collection 
and reject background terms, they do not teach or suggest that the dither and the differential 
data collection must be performed at high-speed. In fact, Leitgeb et al. (including Fercher) 
teach that simply integrating the array detector for a short period of time rejects 1/f noise. 
R. Leitgeb. et al.. Optics Express, 11. 889-894 (21 April 2003). They explicitly state that: 
"Fourier domain systems record one full A-scan In parallel. For short exposure times (<1 ms) 
also In this case 1/f noise will be neglected." This Is Incorrect To first order, a single short 
exposure passes all frequencies from DC up to approximately the reciprocal of the pulse width. 
To reject DC and low frequencies, the array must be integrated for a short period time, and the 
time separation between the shifted interferograms must be short, for it is this time separation 
— not the integration time — that determines the measurement bandwidth. In other words, 
synchronous detection as in the present invention must be used. 

BRIEF SUMMARY OF THE INVENTION 
The present invention is of an apparatus and method for spectral interferometry, 
comprising: providing an interferometer comprising a light source; and dithering an element to 



provide a continuous relative phase shift between target and reference arms of the 
interferometer. In the preferred embodiment, the light source is mode locked. The invention 
subtracts spectra that differ by the phase shift to create a differential spectral interferogram. 
then performs a Fourier transform of the differential spectral interferogram. whereby an 
amplitude of a signal of interest is improved by a factor of approximately two as compared to 
non-differential spectral interferometry. Both real and imaginary components of a complex 
Interferogram are determined, which substantially removes 1/f noise. The element preferably 
comprises a piezo translator. Lock-in detection for each of one or more pixels of a detector is 
preferably employed. The phase shift is preferably approximately tt. the invention substantially 
rejects 1 /f noise, substantially rejects low frequency noise, and substantially reduces detection 
bandwidth. Synchronous detection is preferably employed, most preferably lock-in detection. 

The invention is also of an apparatus and method for spectral interferometry. comprising: 
providing an interferometer comprising a light source; dithering an element to provide a 
continuous relative phase shift between target and reference amis of the interferometer; and 
providing synchronous detection for each of one or more pixels of a detector. 

The invention is further of an apparatus and method for spectral interferometry, 
comprising: providing an interferometer comprising a light source; employing an element to 
provide a relative phase shift between target and reference arms of the interferometer; and 
providing synchronous detection for each of one or more pixels of a detector. 

Objects, advantages and novel features, and further scope of applicability of the present 
invention will be set forth in part in the detailed description to follow, taken in conjunction with 
the accompanying drawnngs. and in part will become apparent to those skilled in the art upon 
examination of the following, or may be learned by practice of the invention. The objects and 
advantages of the invention may be realized and attained by means of the instrumentalities and 
combinations particularty pointed out in the appended claims. 



BRIEF DESCRIPTION OF THE DRAWINGS 
The accompanying drawings, which are incorporated into and form a part of the 
specification, illustrate several embodiments of the present invention and, together with the 
description, serve to explain the principles of the invention. The drawings are only for the 
purpose of illustrating a prefen-ed embodiment of the invention and are not to be construed as 
limiting the invention. In the drawings: 

Fig. 1 is the schematic diagram of the preferred DSI Instrument according to the 
invention. The instrument is preferably constructed from the general form of a Michelson 
interferometer. The sample is in one arm and a dithered mirror is placed In the other ami. The 
output is spectrally resolved. Detection Is synchronous with the dither. 

Fig. 2 is the interferograms (insets) and their Fourier transfonns that produce one- 
dimensional depth profiles of a test sample constructed of microscope glass cover slides: (a) SI 
Image with zero displacement of the reference arm mirror; (b) SI Image with a n phase shift in 
the reference arm; (c) DSI image. The light source is a superluminescent diode centered at 
--840 nm with a bandwidth of 12.3 nm. 

Fig. 3 are the cross-sectional (depth vs. lateral position) images (500x256 pixels) of a 
region of a Xenopus tadpole head obtained by DSI without averaging (a) and with averaging of 
100 Fourier-transformed DSI spectra for each lateral position (b). The light source is the 
superluminescent diode. 

Fig. 4 is the DSI spectra of test samples constructed of glass slides obtained with the 
Ti:sapphire (a) and super-luminescent diode (b) light sources. The FWHM of the peaks are 
indicated near the peaks. 



Fig. 5 is the schematic diagram of a test sample constructed of a RonchI grating and a 
microscope cover slide (a) and the corresponding cross-sectional DSI image (b). The light 
source Is a Tiisapphire laser. 

Fig. 6 is a plot showing the reciprocal of the noise amplitude versus square root of the 
averaging number (points) and a linear fit to the experimental data (line). 

Fig. 7 is a plot evaluating the dynamic range of the DSI imaging instrument. Dashed lines 
indicate the noise levels for averaging of the signal over 1 and 1 00 scans. The light source is a 
Ti:sapphire laser. 

Fig. 8 is a DSI spectrum of a mirror located 240 pm in front of the reference plane 
obtained with true lock-in detection (a) and its Fourier transform (b). The light source Is a 
passively mode-locked Ti:Sapphire laser. The fitted Gaussian peak width is indicated. 

Fig. 9 shows the plot used to determine the dynamic range of the DSI imaging instrument 
with scanning monochromator/lock-in detection. The plot shows the intensity of the peak 
corresponding to the image of the mirror with different neutral density (ND) filters inserted in the 
sample arm of the interferometer. The dashed line indicates the noise level. The light source 
is a Ti:sapphire laser. 

DESCRIPTION OF THE PREFERRED EMBQDIMENT.q 
(BEST MODES F OR CARRYING OUT THE INVENTION) 
The present invention Is of a differential spectral interferometry (DSI) method and 
apparatus which combines the advantages of both OCT and SI. The invention Includes a 
modification of Si that requires only a small m dither (n phase shift) of the interferometer's 
reference path or. equivalently. a m or id2 phase shift in the reference arm. Although DSI 
does not distinguish between the positive and negative optical path differences, it eliminates the 
DC and autocorrelation terms. It Is much simpler than a five-frame complex reconstruction 



method and is not sensitive to the sample Instability. A differential spectral interferogram is 
obtained by subtracting two spectra that differ only by a 7t phase shift in the reference path. A 
Fourier transform of the differential spectral interferogram produces an Image that is free of the 
background inherent In SI, with the signal of interest increased by a factor of two (compared to 
SI). As with SI, the spectral data are collected in a parallel way, however the data acquisition 
can be performed in an AC mode, which is compatible with lock-in detection, allowing low- 
bandwldth measurements. This results in a substantial improvement of the imaging dynamic 
range, making It comparable to OCT, while maintaining the simplicity of SI. 

Consider two plane electromagnetic waves Interfering at the output of the 
interferometer. The resulting spectrum is: 

S(a}) =1 E^ia?) f + 1 62(0?) |2 +2E^(a})E2{o})cos(^^(co) - </>2{co) - cor) (1) 

where ^yis the frequency, S{co) is the spectrum at the output of the interferometer, E^(eD) and 
E2(q}) are the magnitudes of the electric fields in each arm of the interferometer, 4>^{a^ and 
are the phases of the two waves, and r Is the optical delay between the two arms of the 
interferometer. The first two terms are the DC spectral terms from the light source. The third 
term contains the desired information. 

If one adds a small time delay, Ar, and subtract that spectrum from the spectrum in 
Eq. (1), we get the following difference spectrum SdiKflj): 

Sdif (fi^) = ^iM^2(fi^)[cos{^i(iy) - ^2(^) ~ ^r) - COS{^^{a)) - ^2(0^) - fi?(r + A r))] (2) 

By applying familiar trigonometric identities, Eq. (2) can be simplified to the following form: 



^dif (^) = ^^{o))^2i^){^^^{4>^{co) - (p2{o>) - (or) - oos{4>^{q}) - - fi;r)cos(ffl!Ar) 

- Sln(^i(6?) - ^(CO) ~ 6}T)s\X\{a)^T)] 



If Aris set such that oAt^ n (Ar is reasonably assumed to be much smaller than the inverse of 
the light source bandwidth), then Eq. (3) reduces to: 

Following from the above description, when a n phase shift modulation is introduced 
into the reference arm. the DC and autocorrelation terms are cancelled and there is two-fold 
increase of the useful signal amplitude. 

A schematic diagram of an apparatus 10 according to the invention used to obtain 
images of test targets 22. such as biological samples, is shown in Fig. 1. The DSI instrument is 
based on a Michelson-type interferometer. A laser 12 producing a broadband output, 
preferably either a mode-locked Tirsapphire laser (A = 800 nm, FWHM = 27.3 nm) or a 
superluminescent diode (A, = 840 nm. FWHM = 1 2.3 nm) is used as the light source, which light 
passes through mode matching optics 24 to beamsplitter 26. The interferogram is spectrally 
resolved 14. such as with a 1/4 meter imaging spectrograph, and detected via detector 16, 
preferably with an array detector such as a 1 6-bit CCD camera. The mirror 1 8 in the reference 
arm of the interferometer is dithered, preferably with a piezo-translator at /= 70 Hz. to introduce 
a 71 phase shift to the reference beam. To get a single depth profile, the spectra at the two 
extreme positions of the dithered mirror are acquired, and a Fourier transform is performed on 
the difference spectrum via computer 20. While this detection system does not completely 
eliminate the dynamic range problem inherent in CCD's, it demonstrates the basic utility of DSI. 
To fully test DSI. lock-in detection experiments were also performed (see below). 

Figure 2 compares the one-dimensional SI scans (obtained with the phase shifts of 0 
and 7t) to the DSI scan acquired for a test sample constructed of microscope glass cover slides. 
It is seen from Fig. 2 that DSI with differential detection provides a very efficient suppression of 
the background features, including the DC and autocon-elation terms in the vicinity of the 



position of the reference plane (zero depth). As expected, the intensities of the DSI signal 
peaks increase about twice their original value in the SI scans. 



The axial resolution is 1 1 with the laser and 24 Mm with the SLD light source; the 
lateral resolution Is 10 urn. As in SI, the axial resolution is limited by the bandwidth of the light 
source, while the lateral resolution is determined by the focusing conditions (f = 2.52 cm, in 
these experiments). The free-space depth range is limited by the spectral resolution of the 
spectrograph, and is about 0.3 cm for the present instrument. The dynamic range, measured 
by inserting different neutral density filters in the sample ami of the interferometer, was found to 
be about 75 dB. Unlike SI. DSI removes the background features allowing the image signal-to- 
noise ratio, which is determined by the random noise of the baseline, to be improved by 
averaging multiple images. By averaging 100 differential spectra, a dynamic range of 94 dB 
was obtained. Further gains in dynamic range result from optimization of the detection system. 
Figure 3 illustrates the improvement in the quality of a DSI image (500x256 pixels) of a region 
of the head olaXenopus Laevis (African clawed frog) tadpole by averaging over 100 Fourier- 
transformed differential spectra for each scan. With the averaging, the effective imaging rate 
was 5.6 ms/pixel. 

To demonstrate the advantages of using lock-in detection, experiments were 
performed with a scanning monochromator and single-channel lock-in detection with a dithering 
frequency of 1.4 kHz. A significant enhancement in the dynamic range up to 105 dB was 
observed. With the use of a photodiode array detector with true lock-in detection for each 
individual photodiode. improvement of the image acquisition rate to «10 Hz for video-rate 
imaging is possible. For example, to achieve a 10 Hz rate for a 2D-image consisting of 100 
ID-depth profiles, the acquisition rate should be one ID-scan per 1 ms. There are no 
fundamental limitations in achieving that rate. At present, with one-channel lock-in detection 
we obtain high-quality differential spectra with a dithering frequency of up to approximately 5 
kHz and time constant of the lock-in amplifier of approximately 300 us. Note, for an imaging 
rate of one ID-scan per 1 ms, the noise equivalent bandwidth will be «1 MHz for OCT and «1 



kHz for DSI. This means that, theoretically, with the same power of the light source, quantum 
efficiency and dynamic range of the detectors, multi-channel DSI can provide a signal-to-noise 
ratio (SNR) of about a factor of 30 higher than OCT. assuming the noise scalar as the square 
root of the bandwidth. 

The invention provides a simple and powerful new modification to SI that provides 
effective background suppression and greatly improves the obtainable dynamic range. For 
biomedical applications, the multi-channel version of DSI is competitive with OCT in image 
quality and acquisition rate, while being more robust and mechanically simpler. 

Differential spectral interferometry preserves the features of OCT and SI alone while 
providing the following advantages: 

DSI is optically and mechanically simple and reliable. DSI is based on a fixed 
interferometer - the reference arm requires only a small dither of a few microns rather than 
several millimeters greatly simplifying alignment and design. 

DSI is inherently parallel leading to faster imaging rates. DSI uses frequency domain 
interferometry allowing the entire interferogram to be recorded simultaneously. 

DSI reduces system bandwidth requirements resulting in clearer images. The inherent 
parallel data acquisition of DSI reduces the bandwidth of the detection electronics from a few 
MHz to a few 10's of Hz. 



DSI hardware is similar to OCT. speeding regulatory approval. 



In determining the utility of tlie invention, the following questions present themselves: 

(1) Can DSI provide sufficient spatial resolution and signal-to-noise ratios for imaging 
tissue microstructures? Specifically, are axial resolutions near 10 ^m achievable by this 
method? 



(2) Can DSI provide sufficient contrast for cross-sectional tissue imaging? Specifically, 
can a high resolution CCD camera (or other array detector) provide the necessary dynamic 
range for imaging to depths of a few mm in scattering media? 

(3) How does imaging DSI compare with current coherence imaging methods? Will video 
rate imaging ultimately be achievable in a commercial instrument? 

In testing, questions (1) and (2) were answered affirmatively, i.e., It was demonstrated 
that differential spectral interferometry can be used to image biological tissues. In addressing 
question (3), DSI can be used for high-speed acquisition of cross-sectional images of biological 
objects with the required resolution and dynamic range. The image quality is comparable to 
that of OCT, and there are no fundamental limitations which will prevent acquiring images at 
video rates. 



Optical Coherence Tomography 

Optical coherence tomography (OCT) is the application of optical coherence 
reflectometry (OCR) to cross-sectional Imaging in biological tissues. Optical coherence 
reflectometry is a technique that has been used for many years for high resolution phase and 
distance measurements, and surface contour mapping in such materials as optical 
components, eye structures, and integrated circuits. Since 1991. development and 
demonstration of OCT for a variety biomedical applications has been pursued. 
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Optical coherence reflectometry (or tomography) is essentially a time domain technique 
where distance is measured using optical delay times. A spectrally broad, low coherence light 
beam, such as output from a femtosecond laser or super luminescent diode (SLD). is split into 
sample and reference beams within a Michelson interferometer. The sample beam is directed 
into the specimen, while the reference beam travels a path of nearly equal distance. The 
reference arm length is rapidly modulated. Light retuming from the sample recombines with 
the reference beam at the beam splitter. A single element. AC detector is used to detect 
interferometer output. Since the light has a short coherence length, interference occurs only 
when the optical time-of-travel through each arm of the Interferometer is the same. The 
interferometer pattern yields the spatial profile of the sample directly. 

Optical coherence interferometry techniques, such as OCR and OCT. actually measure 
differences in phase between the reference and sample amis. Different optical path lengths 
between the two arms, whether resulting from an actual difference in distance traveled, or a 
difference in index of refraction, result in a phase shift between the two beams. This phase 
shift is recorded in an interference pattern. Knowing the length of the reference arm and the 
index of refraction of the sample allows determination of the depth from which light is reflected 
or backscattered in the sample. Rapidly scanning the reference arm length allows signal from 
varying depths in the sample to be collected, while scanning the beam across the specimen 
provides imaging in the transverse plane. Multiple axial scans are made at a series of lateral 
positions to provide a 2-dimensional cross-sectional map of reflection or bacl<-scattering 
locations In the sample. 



Most high resolution optical imaging techniques, such as confocal, two-photon and 
near-field methods, require high intensity laser sources. In confocal and multi-photon 
microscopes, the longitudinal (depth) resolution is detennined by the beam focus and collection 
apertures. For OCT the longitudinal resolution is inherent in the light source - it is determined 
by the coherence length. This means that a super luminescent diode (SLD) can be used while 
still maintaining good resolution; 15 urn depth resolution has been demonstrated with an SLD 



source. As with other optical imaging methods, transverse resolution is determined by beam 
focusing and apertures. 



OCT Is often implemented with optical fibers for both interferometer arms. This mal<es 
the technique suitable for in vivo and endoscope applications. Other advantages are that low 
powers (micro- to milliwatts) and near infrared wavelengths, where tissues and cells are most 
transparent, are used. This avoids photo-induced damage to the specimen. In the past few 
years, research into OCT applications including endoscopic OCT instalments, Doppler OCT for 
blood flow imaging, in vivo measurements in several different organs, and work towards 
disease and cancer detection have been reported. Other developing applications include 
assessment of tissue during microsurgery, and imaging of embryonic and neural tissue. 

One drawback to OCT is that real-time imaging is not simple to achieve. The need for 
a high-bandwidth (up to a few MHz) detector precludes the use of cameras at high frame rates. 
Fast imaging also requires high speed modulation of the reference arm, which is not feasible 
using inexpensive mechanical components. New advances such as a Fourier domain optical 
delay line have made video rate OCT possible, but high bandwidth electronics are still required 
even when parallel detection is used. DSI is a frequency domain technique, making it 
inherently parallel. As a result, the bandwidth of the electronics is greatly reduced, lowering 
noise. 



Spectral Interferometrv ^ 

Spectral interferometry (SI) has been applied to high resolution distance and phase 
measurements in a variety of applications. However, until recently SI has received little 
attention in the biomedical optics field. One notable exception is a 1998 publication by HSusler 
and Lindner on dermatological applications of SI, which they term "Spectral Radar." G. 
Hausler. et al.. "Coherence radar and spectral radar - new tools for dermatological diagnosis", 
J. Biomed. Opt. 3. 21-31 (1998). Other biomedical applications have been in ophthalmology. 
A.F. Fercher. etal., "Measurement of intraocular distances by backscattering spectral 
interferometry", Opt. Commun. 117. 43-48 (1995). 
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Like OCT. SI is an optical technique that uses an interferometer and a low coherence, 
spectrally broad light source. However, with SI. the interferometer is fixed; instead of scanning 
the reference arm. the interferometer output is spectrally dispersed before detection. The 
signal from the sample consists of a superposition of monochromatic waves returning from 
different depths in the sample. These interfere with corresponding components from the fixed 
reference arm creating a spectral interferogram. The interferometer output is sent through a 
spectrograph to resolve the interference fringes. 

Depth information is encoded in the frequency of the fringes and can be extracted 
using simple Fourier transform that can be performed in only a few microseconds. With SI. all 
depth information is acquired simultaneously in an inherently parallel manner using a CCD 
array, eliminating the need for axial reference arm scans. The spectrum of the interferometer 
output (sum spectrum), Isi(u)). is given by 



IsM = Ir(^) + /s(<y) + 2^Ir{<o)Is{<») X cos{q>s(to) - <pfi(co) - cot) 



(5) 



where lsi((o) is the interferometer output spectrum, expressed as intensity as a function of 
optical frequency w. Ir((o) is the spectrum of the light returning from the reference arm while 
ls{co) is the spectrum of the light returning from the sample arm. (Pr(co) is the reference arm 
phase, <ps(<o) is the sample arm phase, and t is a fixed optical delay between the two arms. 

The first two terms in Eq. 5 are simply the reference and sample arm spectra. The 
third term - the interference between the reference and sample arms - contains the depth 
(phase) information. As with OCT. a phase delay can result from both differences in distance 
between the reference arm and backscattering loci in the sample and from differences in index 
of refraction. It is necessary to know the sample index of refraction in order to quantitatively 
determine depth. 



10 
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To facilitate extracting the depth information from the SI output, the delay, x. is chosen 
to yield fringes in the sum spectrum, i.e.. so that the optical paths are nearly, but not exactly, 
equal. This delay then remains fixed, and optical path differences resulting from back- 
scattering from various depths in the sample are encoded in the interference fringes. The 
spectral fringes have a period inversely proportional to the optical path difference between the 
two beams. 

The first step is to subtract the individual spectra. Ir((o) and ls((o). from the total signal in 
order to isolate the spectral interferogram. S(co), where 



S(0)) = 2yjlff{o>)ls{0)) X COS(«»s(<») - q)ff{0}) - tDT) (6) 

By Fourier transfonning the result, S(co). one obtains 

3"'[SM] = ^(r-r) + f(-/-r)(7) 



where f(t) is the correlation product between the reference and sample fields. In this example, 
reflection from a single depth in the sample results in a simple function of intensity versus time 
delay. Note that "mirror images" off(t) occur at positive and negative delays. Ifthe time delay 
results purely from reflection or scattering, then time corresponds directly to distance. In tissue 
samples, the index of refraction must also be taken into account. 



Is(<a) « Ir(os) when there is no significant absorption that changes the spectral profile of 
the light returning from the sample (except for possible attenuation at all frequencies due to 
25 multiply scattering media). In this case, subtraction of Is(co) and IrC©) before the Fourier 
transform is feasible. However this is not necessary. If the individual spectra are left in the 
interferogram. they will appear as a background feature centered around zero in the time 
domain and generally can be removed by appropriate filtering. However, the dynamic range 
issue still remains a problem. 
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Complex, multi-component signals are also observed for biological tissues. In highly 
scattering samples, it may be necessary to use a weighting function (typically exponential) to 
account for attenuation due to multiple scattering of light returning from greater depths. 
Otherwise, the mathematical treatment remains relatively simple, requiring only Fourier 
transformation of the spectral interferogram. 

Different ial Spectral Interferometry 

DSI is a modification and improvement of SI that requires only a small dithering of the 
reference path of the interferometer. A differential spectral interferogram is obtained by 
subtracting the individual spectra acquired at two slightly different positions of the dithered 
mirror. Instead of a fixed mirror in one arm of the interferometer, a mirror is used that can be 
dithered slightly (a few microns at most). To complete the subtraction at the dither frequency, 
synchronous detection (preferably lock-in detection) is used, eliminating the spectral 
contributions from the individual pulses that improves the dynamic range. A Fourier transform 
is performed on the differential spectrum to obtain the sample depth profile. By optimizing the 
dither to a TT phase shift, the spectral profile of interest coherently adds while the deleterious 
background is efficiently removed. 

Consider two plane electromagnetic waves interfering at the output of the 
interferometer. The spectrum of the two electric fields is given by: 

S(a,) = |E,(c,)e'<*'<->-'-'i + C.C. + E2(«)e"^''")-''^= + ccf (8) 

where S(u)) is the spectrum. E(w) is the magnitude of the electric field (and hence real), cp(aj) is 
the phase as a function of w. and t is the time and c.c. is the complex conjugate. Completing 
the square yields: 



S(w) = \E^(a>)f + 1£2(<»)|^ + Ei(o)E2(a))(e''*''"'-''^<'">-''^'i-'2) + giM<»)-i^M-ia,it2-t,)^ 



(9) 
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which gives: 

S{o?) = \E,(cof + lE^icof + E,{co)E2{a,)cos(Mco) - ^(o,) - a^r) (10) 
where t is the difference between t, and tz- 

If one adds a small time delay. At. to the time delay between the electric fields, x. and 
subtracts that spectrum from the spectrum in Eq.10. one gets: 

E^(a})E2(a})[cos(MQ}) - ^^(o)) > mr) - OOS{(P^(co) - cf^^^co) - co{r + Ar)](1 1) 



Using: 



cos(<9i - ^2) = cos cos - sin^^ sin 6^2 (12) 



one gets: 



CO%{<t>^ {CO) ~ <j>2 {oj) ~ or) - cos(^i (o)) - ^2(^) - 0}t) COS(fi;A rj 
- s\n(i^^(o)) - ^2(o>) - eoT)s\n{coAT) 



(13) 



If At is set such that wAx - tt. then the above equation reduces to: 



2E^(co)E2(o}) cos{^^{o}) - ^2(^) ~ o^r) (14) 



Thus, the signal of interest. Ei(u))E2(u)). coherently adds while the background terms 
subtracted. The experimental work demonstrates this. 



are 
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Even when wAx is close to tt. equation 13 can be approximated as: 

E,(o)E2(co)[2cos{Mo>) -<^{a>) - eor) - s\n{^{a>) - ^{a>) - a,r)sin(a«^r)](15) 

The second term in equation 8 acts a modulation on the first term producing side bands on x 
spaced by Ax. As long as the resolution » A/2, the spacing of the side bands cannot be 
resolved and does not affect the resolution. For very wide bandwidth sources ~iOo/1 0 or more 
(where Wo is the center frequency), equation 13 cannot be approximated. 

When a tt phase shift modulation is introduced into the reference arm. the DC and 
autocorrelation temis (included in the |E,(u,)|^ and |E,(a,)|^ temis in Eq. 1 0) are cancelled, and 
there is two-fold increase of the useful signal amplitude. The differential spectrum can be 
obtained in the AC mode, which greatly reduces the noise equivalent bandwidth. With the use 
of a photodiode array detector with lock-in detection for each individual photodiode. the 
differential spectra can be collected in an inherently parallel way. Both the reduced noise 
equivalent bandwidth and the multichannel character of data acquisition lead to an increase of 
the signal-to-noise ratio. Our estimates show that multichannel DSI may provide a signal-to- 
noise ratio of a factor of 30 higher than OCT. In addition, if the full dynamic range of the 
detectors can be used, an additional improvement in the signal-to-noise ratio of N"^. where N is 
the number of detectors, may be obtained. 

Taking the difference signal between a small dither in the relative phase between the 
arms of the interferometer is not the only way to reject the DC background, autocorrelation 
terms and 1/f noise. According to J. Schmit. et al., "Extended averaging technique for 
derivation of error-compensating algorithms in phase shifting interferometry." /»pp/. Opt. 34. 
3610 (1995). arithmetically combining several different interferograms with several different 
phase differences can be used to detemnine either the real (cosine) part of the imaginary (sine) 
part of the interferogram. which vwll reject DC background and autocorrelation temis. However. 
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it is also important to measure the component interferograms (that are arithmetically combined 
to obtain the real and the imaginary parts of the interferogram of interest) in a manner that 
rejects 1/f noise. Thus, the component interferograms must be obtained as close to 
simultaneously as possible. 

For example, suppose the interferogram of interest is the arithmetic combination of six 
different interferograms. each labeled I, through le and each with its own relative phase, and 
the interferogram of interest is given by -I, + 51^ +213-101, +3ls -He. The lowfrequency cutoff is 
determined by the amount of time required to measure the set of six interferograms. If 
averaging is required, complete sets of the interferograms should be averaged, which 
increases the low frequency cutoff (the order in which the component interferograms is taken is 
not important) so that it has a high enough frequency to reject 1/f noise. Averaging each 
component interferogram separately lowers the frequency of the low frequency cutoff (makes 
the low frequency cutoff closer to DC) and. therefore, does not reject the 1/f noise. 

Consequently, taking the data in the order l„ 1^, I3, |,. i^. ,3. ,^ ,^ ^ 

Is. Is is not the same as taking the data in the order I,. U, l„ |,. |„ I3. ,3. ,3. 13. ^_ ,^ 
Is. le. In the first case, the low frequency cutoff is approximately given by the reciprocal of the 
time required to measure the first six interferograms. which we will call t. The low frequency 
cutoff is then approximately 1/t. In the second case, the low frequency cutoff is approximately 
the reciprocal of the time required to measure all 1 8 interferograms. which is St. Thus, the low 
frequency cutoff is a factor of 3 smaller in the second case. If the first case is used to average 
the interferogram. all the frequencies between 1/t and DC are rejected. However. In the second 
case, only the frequencies between 1/(3t) and DC are rejected-all the noise in the frequencies 
between 1/(3t) and 1/t is integrated into the measurement with no gain in signal. 



Note that the present invention can also be used with the variation of spectral 
interferometry in which the broadband light source is replaced with a swept wavelength source. 
Because the wavelength is swept, the need for a spectrometer is removed and a single 
30 detector can be used instead of an array detector, which simplifies the device. The 
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disadvantage Is that the data collection time is increased because the data for each wavelength 
cannot be taken in parallel. 



Imaaina spectral in terferometer 

A spectral interferometer was constructed using a Melles Griot 03 BSC 029 cube 
beamsplitter (BS). two Newport 1 0D520ER.2 silver-coated pyrex flat mirrors, and two Newport 
GPX043LR. 16 GRADIUM singlet lenses (f = 2.52 cm). The spectra were analyzed using a 1/4 
meter Chromex imaging spectrograph combined with a cooled 16-bit CCD camera (C7041. 
Hamamatsu). Spectral interferometry was performed using a broad band, low coherence 
length, passively mode-locked Tirsapphire laser (KM Labs) ora super luminescent diode (LDN- 
16. Volga Technology). The emission spectra of the laser and the super luminescent diode 
(SLD) have bandwidth FWHM of I.S.IO" Hz (27.3 nm) and 0.6xl0» Hz (12.3 nm). 
respectively. The instrument was constructed to use both types of light source. The schematic 
diagram of the instrument is shown in Fig. 1. 

The mirror in the reference arm of the interferometer ("Dithered mirror" in Fig. 1) is 
preferably mounted on a ThorLabs KC1-P2 piezo translator, which allows high-precision 
translation of the mirror within up to 5 nm along the optical axis using a ThorLabs MDT693 
piezo controller. 

The sample holder ("Test target" in Fig. 1 ) is mounted on a 9065-XYZ New Focus XYZ 
(3 axis) translation stage. The position of the sample in the XY-plane (perpendicular to the 
laser beam) is controlled by two Newport CMA-12CCCL motors. The CCD-camera, motorized 
XYZ translation stage, and the controller of the piezo translator are interfaced to a computer 
using software created in National Instruments LabView® graphical development environment. 



To evaluate the spatial resolution of the instrument and to obtain images of test 
samples and biological tissues, we used high-speed differential detection. The CCD camera is 
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capabteofacquiring spectra ataspeedof.pto160 Hz in binned mode. In these experiments 
the CCD camera readout rate was set at 140 Hz. Thus, the delay line was dithered at 70 Hz 
synchronized to the readout of the CCD camera. To get a singie depth profile, the spectra at 

the two extreme positions of the mirror were acquired, and the Fourier transform was 
5 performed on the difference spectrum. While this detection system does not compietely 

eliminate the dynamic range problem inherent in CCDs, it demonstrates the utility of DSI. A 

true lock-in detection scheme is discussed below. 

Figure 4 shows an example of experimental SI and DSI spectra and the resulting depth 
profilesofasampleconstructedofmicroscopecoverslides. Thetwo SI spectra correspond to 
a TT phase sm introduced to the reference beam. Both 81 profiles look similar, shovwng the 
peaks at 400 pm. 550 pm. 715 pm, and 1 120 pm that correspond to reflective surfaces of the 
sample. Besides the useful signal peaks, the 81 p^fiies contain strong background, which 
gains oven^helming amplitude at depth close to zero (DC tem,). As seen from Fig. 4. DSI 
provides a very efficient suppression of the background (including the DC temi). Moreover, the 
intensities of the DSI signal peaks increase about twice their original value in the SI scans 
These features of DSI. which could be anticipated from the theoretical description of the 
method (see previous section on DSI), provide dramatic improvement of the SI technique in 
terms of the signal-to-noise ratio and image quality. 

Spatial resolution 

The axial and lateral resolution of DSI images are basically the same as those of 81 
Namely, the axial resolution is determined by the wavelength and bandwidth of the light source, 
and the lateral resolution is detemiined by the incident light focusing conditions. 

Figure 4 shows representative examples of depth profiles of test samples constructed 
ofgIassslides.obtainedwiththelaser(4a)and8LDIightsource(4b). The peaks were fitted by 
the Gaussian function. The FWHM of the fitted peaks are indicated in the figures. The FWHM 
of the peaks were found to be 20 ± 2 pm and 24 ± 2 pm for the laser and SLD light sources, 
respectively. Theoretical bandv^dth-llmited peak widths are »10 pm for the laser (bandwidth 
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1.3x10" HZ) and .22 for the SLD (0.6.10" Hz). For the SLD light source the 
experimental peaK width is Cose to the band«.dth-lin,ited value. Theoretical resolution of the 

.aser light source and theobtainedresolutionare different because the spectral,«.ndowof the 
spectrometer does not allow detection of the whole emission bandvwdth of the laser The 

5 spectral «.ndow is detem,ined by the spectrograph/CCD camera spectrum analyzer, whose 
perfom^ance is optimized to achieve the best spectra, resolution. High spectra, resolution 
provides good imaging depth (a.most 3 mm in free space, v^th this instrument) atthe expense 
Of spatia. rescution. The use of a different (.ower resolution) spectrograph would improve the 
axial spatia. resCution. however the imaging depth would be worse. As seen below, with true 

• lock-in detection, the bandv«dth-.imited axia. resolution is achievable with the Ti:sapphire laser 
light source. 

Figure 5 shows an example of a cross-sectiona. scan of a test samp.e constructed of a 
Ronchi grating (40 .ines/mm) and a microscope cover g.ass slide (thickness 1 10 pm) The 
glass Slide was attached to the grating through spacers (thickness 80 pm). The sample is 
Shown schematically in Fig. 5(a). The DS, 2D-scan produces a clear image of the surface of 
the Ronchi grating (at 424 pm) and two surfaces corresponding to the glass slide (at 176 pm 
and 344 pm). The spacing between the glass slide and the grating, as well the thickness of the 
glass Slide, evaluated from the DS. image, agree well with the corresponding thicknesses 
measured independently «.th a digital caliper (note, the thickness of the glass slide obtained 
from the image. 168 pm. is the actual thickness multiplied by the index of refraction of the 
glass, n «1.5). 

Evaluation of the dy nan nic range 

To evaluate signal-to-noise (S/N) ratio of the invention, a method similar to that 
described by Hausler and Lindner was employed. One detem^ines the minimum intensity of 
l.ght reflected from the sample (relative to the intensity of the incident light) that can be 
Observed v^th the S/N ratio equal to unity. In other words, the dynamic range of the setup is 
evaluated. 
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First, consider the noise. As mentioned above, DSI efficiently suppresses the 
background that is inherent to SI imaging. With SI, the imaging dynamic range is limited by the 
background rather than by noise, which prevents improvement of the Image quality by 
5 averaging multiple scans. To check if the random noise limit can be achieved in DSI, the noise 
of the baseline of DSI images of a mirror was measured v^th averaging over 1 to 100 scans. 
The results are shown in Fig. 6. It is seen from Fig. 6 that the reciprocal of the noise scales 
linearly with the square root of the number of averaged scans. This indicates that the DSI 
dynamic range is limited by random noise, and that the quality of the image can be efficiently 
1 0 enhanced by averaging. 

To measure the dynamic range, an aluminum-coated mirror was used as a sample. 
The intensity of the peak corresponding to the image of the mirror was measured with different 
neutral density (ND) filters inserted in the sample arm of the interferometer. Figure 7 shows the 

1 5 peak intensity plotted versus the ND filter optical density (OD). Because the light in the sample 
arm of the interferometer passes through the filter twice, the OD was calculated as twice the 
optical density of each ND filter. For each filter, OD was determined experimentally by 
measuring the Ti:sapphire laser light attenuation. The dashed lines in Fig. 7 indicate the noise 
level with the signal averaged over 1 and 100 scans. The experimental points in Fig. 7 can be 

20 fitted by a straight line, which intersects the shot noise level at OD = 7.4, thus yielding an 
estimate of the dynamic range of the detection system, equal to 74 dB. This result means that 
one can detect the intensity of light reflected from the sample that is as low as 4x1 0'° of the 
incident light intensity. This is comparable to the dynamic range data for analogous spectral 
interferometry imaging setups.^®'^® Unlike SI. the imaging quality of DSI can be improved by 

25 averaging. As seen in Fig. 7, by averaging over 100 scans, the dynamic range is increased to 
about 94 dB. 

Contrast and resolution of images in scatterino media 

The effect of light-scattering media on the contrast and axial resolution of DSI images 
30 of tissue phantom samples was investigated. In these experiments, both the Ti:sapphire laser 
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and SLD were used as the light sources. The pseudo-lock-in detection scheme described 
above was used. Two tissue phantom samples were constructed by attaching thin glass slides 
(thickness 110 ^im) to the back surface of the glass plate through spacers of different 
thickness, 160 ^m and 360 ^m for the two samples, respectively. DSI images of these 
5 samples, immersed in light-scattering fluid, were taken. Diluted Intralipid (Phospholipid 
stabilized soybean oil 20% fat emulsion) purchased from SIGMA was used as scattering 
medium. As expected, the light-scattering medium significantly decreases the contrast of the 
images. From this experiment, the following effective scattering coefficient; Ms = 34 percent'^ 
cm*^ was obtained. The obtained [Js is in good agreement with the literature data. See, e.g., 
10 H.G. Van Staveren, et al., "Light scattering in intralipld-10% in the wavelength range of 400- 
1 100 nm," AppL Opt. 30, 4507-4514 (1991). Peak widths were also measured as a function of 
Cd. No effect of the light scattering was observed. 

In combination with the dynamic range measured for our instrument, these data 
15 provide an estimate of the maximum imaging depth that can be achieved in imaging real 
human tissues. A 10% Intralipid emulsion is frequently used for modeling the scattering 
properties of typical human tissues. The effective scattering coefficient of 10% Intralipid is 
=340 cm'\ Combined with the measured 74 dB dynamic range, this yields a maximum imaging 
depth of about 500 ^m. This is a rather conservative estimate of the scattering-limited imaging 
20 depth. DSI data averaging will increase the dynamic range and. therefore, improve the imaging 
depth (up to -640 jxm with averaging over 100 scans). However, even with the 500 |im imaging 
depth, the method is useful for certain important medical applications, e.g., diagnostics of skin 
cancer. Further improvements in detection systems, such as arrays of phase sensitive 
detectors, will further increase the imaging depth. 

25 

Figure 3 shows 2D-images of a region of the head of Xenopus Laevis (African clawed 
frog) tadpole obtained without averaging and with averaging over 100 Fourier transformed 
spectra for each one-dimensional scan. The improvement of the image quality with averaging 
is cleariy seen. The images shown In Fig. 3 indicate that DSI with pseudo lock-in detection is a 
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promising technique for in vivo imaging. Compared to OCT, the main drawbacks of this 
pseudo-lock-ln version of DSI are relatively small dynamic range and a low acquisition rate. 
However, these limitations are not fundamental. Switching to multichannel lock-in detection 
dramatically improves both the dynamic range and acquisition rate, as next discussed, making 
5 DSI competitive to OCT for biomedical imaging. 

Using tme lock-in detection to obtain DSI imaoes 

The utility of DSI has been shown above with the use of pseudo-lock-in detection. A 
fast 16-bit CCD Imager was binned to provide a read out at 140 Hz (dither frequency was 70 Hz 
10 and DSI interferograms were obtained at 70 Hz). However, this was still a hIgh-bandwIdth 
measurement. Each pixel was being read at 250 kHz, providing only a small bandwidth 
advantage over OCT. 

The inherent parallel nature of DSI significantly reduces the bandwidth requirements 
15 over OCT. reducing noise and increasing imaging rates. If one employs an array detector with 
true lock-in detection for each pixel, one can reduce the bandwidth to 30 Hz for video frame 
rates (in contrast to the bandwidth of an OCT device, which would be a1 MHz), However, 
because of equipment limitations in the earlier testing, in order to test DSI using true lock-In 
detection, one must give up the Inherent parallel data acquisition by using a single element 
20 detector. A scanning monochromator was employed together with single channel lock-in 
detection to fully test DSI. Multichannel phase sensitive detection is, however, preferred. 

To obtain the DSI images, a sine voltage waveform was applied to the piezo translator 
of the reference arm mirror to dither it at 1.4 kHz. The wave amplitude was optimized to 

25 produce a tt phase shift. The differential interferograms were recorded by using the Chromex 
imaging spectrometer in the scanning mode and lock-In detection of the signal from a 
photodiode coupled to a laboratory-built translmpedance preamplifier. In these experiments, an 
aluminum-coated mirror was used as a sample. Figure 9 shows an example of a differential 
spectral interferogram and its Fourier transform obtained for a mirror sample mounted In the 

30 sample arm 240 pm In front of the reference plane with the laser light source. 
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With the scanning monochromator/lock-in detection scheme, one has no limitation on 
the working spectral window. As seen from Fig. 8(a). the whole laser emission spectrum is 
detected, therefore, nan^ower peaks can be expected in the measured depth profile. In 
Fig. 8(b), the fitted Gaussian peak width is indicated. With lock-in detection, the peak widths 
were found to fall within a range of 1 1 .3 ± 1 .0 pm. which is close to the Ti:sapphlre laser light 
source bandwidth limit of -10 pm. 

In order to choose an appropriate size of an array detector, a simulation study was 
performed to determine how the peak width depends on the number of points in the DSI 
spectrum. The original experimental DSI spectrum (containing 512 points) shown in Fig. 8 was 
truncated from both ends to leave only 256. 128 or 64 spectral points. Fourier transforms 
perfomied on those spectra produced the peaks with FWHM of 1 1 .2 pm, 13.6 pm. 19.4 pm. 
and 99 pm for the spectra containing 512, 256, 128, and 64 points, respectively. This result 
shows that an array detector containing as few as 128 elements would provide reasonably 
good axial resolution. 

The dynamic range of a DSI Instrument with single-channel lock-In detection was also 
measured. To determine the dynamic range, the technique described above was employed. 
As seen in Fig. 9. there is a dramatic improvement in the dynamic range compared to the 
spectrograph/CCD arrangement (compare Fig. 7). According to Fig. 9, with true lock-in 
detection, the dynamic range is at least 105 dB. The dynamic range can be even larger, if a 
more sophisticated (low-noise) signal pre-amplifier is used. 

Because of the use of the single-channel detection scheme, the acquisition rate Is 
relatively low. For example, it takes about a half minute to record a one-dimensional profile 
shown in Fig. 8. With the use of an array detector with true lock-in detection for each pixel, it is 
possible to improve the image acquisition rate to --1 0 Hz for video-rate imaging. For example, 
with a 128 element array detector, to get a 64x100 image at video rate, the acquisition rate 
should be one ID-scan per 1 ms. which is quite reasonable. In the shot-noise limit, for the 
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same effective acquisition rate. DSI with multi-channel lock-in detection would have S/N ratio 
advantage over OCT of about the square root of the number of channels. 

Summary 

A novel modification of the Spectral Interferometry (SI) imaging method, Differential 
Spectral Interferometry (DSI). has been proposed and evaluated. DSI combines the 
advantages of both Optical Coherence Tomography (OCT) and SI. As with SI. the spectral 
data are collected in a parallel way, however the data acquisition is performed in an AC mode, 
allowing low-bandwidth measurements. The results show that DSI provides a dramatic 
improvement over the SI method in terms of the dynamic range and background suppression. 
With an instrument it was shown that DSI provides a free-space imaging depth of about 0.5-3 
mm (depending on the scattering properties of the sample) with spatial resolution of about 1 1 
pm (which is limited by the bandwidth of the light source). It has been demonstrated that, with 
DSI. the dynamic range of at least 1 05 dB is achievable. With the use of an array detector with 
multichannel lock-in detection, DSI is competitive with OCT in terms of the dynamic range and 
image acquisition rate, while maintaining the simplicity of SI. The design of the individually 
addressed array detector combined with lock-in detection results in an instrument that allows 
high-quality video rate imaging useful for, among others, applications in medicine, biology, and 
material science. 

Extensions of DSI 

By taking the difference signal between a small dither In the reference arm of an 
Interferometer, one directly measures the cosine component of the complex interferogram. By 
shifting the dither point (nominally by one-fourth the center wavelength), the sine component of 
the complex Interferogram can be measured. The cosine and sine components of the complex 
interferogram can be thought of as the real and imaginary parts, respectively. Thus, the 
measurement of the intensity and phase of the entire complex interferogram can be obtained 
using two dithers and a single lock-in detection scheme allowing for real-time measurement of 
the complex interferogram. This removes ghosting and double images in the interferogram 
while allowing for a detection scheme no more complex than required for DSI. One dither has 
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the amplitude required to measure the sine and cosine components while the other dither shifts 
between sine and cosine component measurement. Other dither waveforms can be used to 
obtain the complex interferogram as well. However, some dither waveforms are more 
conducive to lock-in detection than others. Indeed, precalculated waveforms could be 
designed to optimize the signal recovery and compensate for mechanical response. 

To further improve 1/f noise rejection, the magnitude of the Fourier transform of each 
spectral difference (or spectral set, as in the case of the recovery of the complex interferogram) 
can be averaged rather than the spectral differences themselves. In the case of the recovery 
of the complex interferogram, the magnitude of the Fourier transform of each measured 
complex interferogram can be averaged. Thus, rather than averaging the differences, or 
equivalents, the real and/or imaginary parts of the interferogram, the Fourier transform of the 
interferogram (complex or real) is calculated, the magnitude is taken, and the Fourier 
magnitudes are averaged. 

In imaging applications, a wide bandwidth light source is typically used. Obtaining the 
sine and cosine components using a difference method is absolutely accurate only for a single 
wavelength of light. However, measuring the sine and cosine components of the complex 
interferogram using the difference method is typically accurate for imaging requirements even 
though wide bandwidth sources are used. Adjustments can be made to the dithering waveform 
to make the measurement more robust for wide bandwidth sources, if required. Also, full 
digitization of the signal is possible if the waveform is not conducive for measurement with lock- 
in detection. A digital signal processor can process the digitized signal to obtain the fully 
corrected interferogram's intensity and phase. 

By measuring the complex interferogram in real-time, not only are the stability 
problems mitigated, but addition and subtraction of the spectra that comprise the complex 
interferogram measurement efficiently removes background artifacts (as well as other artifacts) 
and allows for shot-noise limited detection of the complex interferogram. 
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This method is an extension of the phase retrieval methods outlined by J. Schmidt and 
K. Creath in their >App//edOpf/c$ paper (Vol. 34. No. 19, pp. 3610-3619). One of the methods 
of Schmidt and Creath are applied to spectral interferometry by WojtkowskI, et al., supra, but 
not In the real-time fashion outlined herein. 

Industrial ADDlicabllitv : 

The invention Is further illustrated by the following non-limiting examples. 

EXAMPLE 1 

The invention can be used to develop a device for determining structural defects in 
composite materials. The device can be placed on composite and the internal structure of the 
composite material can be determined. 

EXAMPLE 2 

The Invention can replace optical coherence tomography In any application that OCT 
can be used. 

EXAMPLE 3 

The invention can be used to guide surgical instruments. 
EXAMPLE 4 

If the complete phase of the interferogram can be determined, then the spectral and 
structural information can be extracted at the same time, within the limits of the uncertainty 
relationship. Thus, the Invention can be used for spectral OCT. This Is a great simplification 
over OCT devices and increases the signal-to-nolse ratio over spectral OCT measurements 
because the inventive technique greatly decreases the detection bandwidth required for 
measurement. 



-30- 

The examples can be repeated with similar success by substituting the generically or 
specifically described reactants and/or operating conditions of this invention for those used in 
the preceding examples. 

Although the invention has been described in detail with particular reference to these 
preferred embodiments, other embodiments can achieve the same results. Variations and 
modifications of the present invention will be obvious to those skilled in the art and it is intended 
to cover in the appended claims all such modifications and equivalents. The entire disclosures 
of all references, applications, patents, and publications cited above, and of the corresponding 
application(s), are hereby incorporated by reference. 



